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ABSTRACT 

We have obtained intermediate resolution optical spectra of the black-hole candidate 
Nova Sco 1994 in June 1996, when the source was in an X-ray/optical active state 
(R ^15.05). We measure the radial velocity curve of the secondary star and obtain a 
semi-amplitude of 279±10kms _1 ; a value which is 30 per cent larger than the value 
obtained when the source is in quiescence. Our large value for K% is consistent with 
60^7 per cent of the secondary star's surface being heated; compared to 35 per cent, 
which is what one would expect if only the inner face of the secondary star were 
irradiated. Effects such as irradiation-induced flows on the secondary star may be 
important in explaining the observed large value for Ki- 

Key words: accretion, accretion discs - binaries: close - stars: individual: X-ray 
Nova Sco 1994 (GRO J1655-40) - X-rays: stars. 



1 INTRODUCTION 

The soft X-ray transients, a subclass of the low-mass X-ray 
binaries distinguished by their X-ray outbursts, have proved 
to be an ideal hunting ground for stellar-mass black hole 
candidates (Tanaka & Shibazaki 1996). The system Nova 
Sco 1994 (=GRO J1655-40) is particularly interesting, since 
as well as being a source of superluminal jets (Zhang et al., 
1994; Harmon et al., 1995), its optical brightness and partial 
eclipse features mean that it is one of the few systems that 
has yielded a reliable estimate for the mass of the collapsed 
star. 

Nova Sco 1994 was discovered on July 27 1994 with 
BATSE on board the Compton Gamma Ray Observatory 
(Zhang et al., 1994). It has been studied extensively dur- 
ing the past few years in X-rays and at optical and ra- 
dio wavelengths (Bailyn et al., 1995a and b, Zhang et al., 
1995, van der Hooft et al. 1998). Strong evidence that 
the compact object in Nova Sco 1994 is a black hole was 
presented by Bailyn et al. (1995b) who initially estab- 
lished a spectroscopic period of 2.601 ± 0.027 days, clas- 
sified the secondary as an F2-F6IV type star and sug- 
gested a mass function /(M)=3.16±0.15 Mq. An im- 
proved value of /(M)=3.24±0.09 Mq was presented by 
Orosz & Bailyn (1997) using both quiescent and out- 
burst data, derived from a radial velocity semi-amplitude 



of 228.2±2.2 kms" 1 . Shahbaz et al. (1999) using only 
quiescent data, determined the true radial velocity semi- 
amplitude 7^2=215. 5±2. 4 km s _1 which gives a revised value 
for the mass function of f(M) = 2.73±0.09 Mq. They also 
measured the rotational broadening of the secondary star 
which then gives the binary mass ratio q ~ 0.39 (=M2/Mi, 
where Mi and M2 are the masses of the compact object and 
secondary star respectively). 

The effect of heating of the secondary is to shift the 
'effective centre' of the secondary, weighted by the strength 
of the absorption lines, from the centre of mass of the star. 
One expects that this results in a significant distortion of 
the radial velocity curve and renders a sinusoidal fit clearly 
inadequate, leading to a spuriously high radial velocity semi- 
amplitude. In order to quantify this effect we have deter- 
mined the radial velocity variations of the secondary star 
in Nova Sco 1994, when it was in outburst and compared 
our results with others obtained using data taken when the 
source was in different X-ray states. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Spectroscopy 

Intermediate resolution optical spectra of Nova Sco 1994 
were obtained on 1996 June 20-24 with the 1.54-m Dan- 
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Table 1. Log of observations 

Night # # of spectra seeing (") Orbital phase 



covered 

Night 1 11 2.0 0.722-0.834 

Night 2 5 2.6 0.100-0.225 

Night 3 11 2.3 0.501-0.592 

Night 4 10 1.6 0.883-0.964 

Night 5 10 3.0 0.270-0.353 



ish Telescope at the European Southern Observatory (ESO) 
in Chile using the Danish Faint Object Spectrograph and 
Camera (DFOSC). We used grating #8 which gave a dis- 
persion of 1.26 A per pixel and a wavelength coverage from 
5865-8336A. The Loral 2048x2048 CCD was used, binned 
by a factor two in the spatial direction in order to reduce 
the readout noise, but not binned in the dispersion direc- 
tion. The seeing during the observations was poor and vari- 
able (see section 2.2) so we used a slit width of 2'.'5 on the 
first night and then 2'.'0 for the other nights. This resulted 
in spectral resolutions of 7.6 A and 5.5 A for the first and 
other nights respectively. Wavelength calibration was per- 
formed using a Cu-Ar arc. A total of 47 spectra were taken 
each having exposure times of 1800s (see Table 1 for details). 

The data reduction and analysis was performed us- 
ing the Starlink FIGARO package, the PAMELA routines of 
K. Home and the molly package of T. R. Marsh. Removal of 
the individual bias signal was achieved through subtraction 
of a median bias frame. Small scale pixel-to-pixel sensitiv- 
ity variations were removed with a flat-field frame prepared 
from observations of a tungsten lamp. One-dimensional 
spectra were extracted using the optimal-extraction algo- 
rithm of Home (1986), and calibration of the wavelength 
scale was achieved using 5th order polynomial fits which 
gave an rms scatter of 0.03 A. The stability of the final cali- 
bration was verified with the OH sky line at 6300. 3A whose 
position was accurate to within 0.1 A. 

2.2 Photometry 

Using the same setup as for the spectroscopy, we also ob- 
tained limited Bessell r-band images of Nova Sco 1994 every 
night. The data were debiased using a median bias frame, 
but not flat-fielded, as none were taken. These images were 
used to estimate the seeing each night (see Table 1). We 
applied aperture photometry to Nova Sco 1994 and several 
nearby comparison stars within the field of view. Johnson 
V- and ii-band magnitudes of these comparison stars were 
made available to us by J. Orosz. We determined the rela- 
tive magnitude of Nova Sco 1994 with respect to three stars 
having a range of colours [(V - i?)=0.49, 0.77 and 1.23]. 
Assuming that the colour correction between the two filter 
systems is small (<0.05 mags; similar to the accuracy of 
our photometry) and that Nova Sco 1994 has a colour in 
the same range as the comparison stars used, we estimate 
R ~15.05 for Nova Sco 1994. 

3 THE SPECTRA OF NOVA SCO 1994 

In Figure 1 we show the variance-weighted average and also 
the nightly averages of the Nova Sco 1994 spectra. A strong 



Table 2. Ha equivalent width 



Night # 


Orbital phase 


EW (A) 


Night 1 


0.78 


9.8±0.5 


Night 2 


0.16 


11.9±0.3 


Night 3 


0.54 


10.6±0.4 


Night 4 


0.93 


3.8±010 


Night 5 


0.32 


10.4±0.8 



Ha emission line (mean equivalent width of 7.5±0.06A) and 
a much weaker Hel 6678A (equivalent width of 0.4±0.03A) 
can be seen. In Table 2 we list the Ha equivalent width for 
the nightly averages. The Fel absorption blend at 6485, 6496, 
6499, 6502A is also visible. These features are used to deter- 
mine the radial velocity of the secondary star (see section 
4). The 6613A diffuse interstellar band is also present. 

The emission lines in Nova Sco 1994 are double-peaked, 
which is presumably a consequence of the system being at 
high inclination. We can compare the observed peak-to-peak 
half separation of the Ha emission line (which arises from 
the accretion disc) with the projected velocity of the outer 
disc edge. In a binary system with a mass ratio > 0.25 it 
is generally assumed that the accretion disc cannot grow 
larger than the tidal truncation radius, (Paczynski 1977; 
Whitehurst 1988; Osaki, Hirose & Ichikawa 1993), which is 
approximately given by r d = 0.60a/(l + q) for 0.03 < q < 1, 
where a is the binary separation (Warner 1995). Given the 
system parameters (P or f,=2. 62168 days; q ~ 0.39; i ~ 69 de- 
grees; Mi ~6.7 Mq see Shahbaz et al., 1999) the minimum 
value for the projected velocity of the accretion disc rim 
is ~394 kms^ 1 . The observed peak-to-peak half separation 
of the Ha emission line (see Figure 1) in late June 1996 is 
385±8 km s -1 (measured by fitting the profile with a double 
Gaussian), which implies that the accretion disc is close to 
its maximum possible size. Soria et al. (1998) estimate the 
Ha half peak-to-peak separation to be <350/2 kms _1 and 
<550/2 kms _1 for their August/September 1994 and June 
1996 Ha observations respectively, velocities much lower 
than expected, suggesting that the Ha emission line arises 
from non-Keplerian regions/flows in the accretion disc. 



4 THE RADIAL VELOCITY OF THE 
SECONDARY STAR 

The radial velocities of the F-type secondary star in Nova 
Sco 1994 were measured from the spectra by the method of 
cross-correlation (Tonry & Davis 1979) with a template star. 
Prior to cross-correlation the spectra were interpolated onto 
a logarithmic wavelength scale (pixel size 55 kms -1 ) using 
a sin x/x interpolation scheme to minimize data smoothing 
(Stover et al. 1980), and then normalised. The template star 
spectrum (HR2906; F6v) was then artificially broadened by 
90 kms _1 (Shahbaz et al., 1999) to account for the rotational 
velocity of the secondary star. Note that the orbital smearing 
of the Nova Sco 1994 spectra through the 1800s exposure is 
at most only 10 kms -1 , much less than the resolution of the 
data. Only regions of the spectrum devoid of emission lines 
(6400-6520A) were used in the cross-correlation. The radial 
velocity of the template star (derived using the position of 
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the Ha absorption line to be — 7 kms -1 ) was then added to 
the radial velocities of Nova Sco 1994. 

Using the orbital ephemeris given by van der Hooft et 
al (1998) we phase-folded and binned the heliocentric radial 
velocities (see Figure 2). From figure 2, it can be seen that 
the radial velocity measurement at phase 0.2 does not fit the 
general pattern of the sinusoidal modulation present in the 
data. This data point was the total of three radial velocity 
measurements taken on the second night (21st June 1996). 
Although the seeing and quality of the spectra taken during 
this night were not as good as the others, no obvious reason 
could be found as to why these spectra gave much lower 
radial velocities than expected. A sine wave fit to all the 
data points does not give an adequate fit (x„=6.9). However, 
removing the discrepant data point and then performing a 
sine wave fit yields a xt °f 1-5, a semi-amplitude K2 = 
279 ± 10 kms -1 , systemic velocity 7 = —155 ± 7 and a 
phase shift of —0.043 ± 0.0050 (l-cr errors are given). We 
also fitted the radial velocity curve with an eccentric orbit, 
but found the fit to be less than 50 percent significant. 



5 THE EFFECT OF IRRADIATION ON THE 
SECONDARY STAR'S RADIAL VELOCITY 

Three absorption line radial velocity curves have been ob- 
tained for Nova Sco 1994, using the same absorption fea- 
tures of the F6iv secondary star and the standard method 
of cross- correlation. However, in each case the system was 
observed to be in a different X-ray state. A sinusoidal fit 
to the outburst data taken in April/May 1995 of Orosz 

6 Bailyn (1998) gives a radial velocity semi-amplitude of 
K oba =230±2 kms -1 . During this period BATSE did not 
detect the source, so we can only put an upper limit of 
2.4xl0 36 ergs -1 (<0.03 photons erg cm -2 s -1 in the BATSE 
20-350 keV energy range) to the X-ray luminosity of the 
source. This upper limit alone does not allow us to state 
unequivocally that the source was not active at X-ray en- 
ergies, but optical observations suggest that the source was 
not in quiescence (V—16.5; Orosz & Bailyn 1998). In sec- 
tion 4 we determined #2=279=1=10 kms _1 from data taken in 
June 1996 when RXTE ASM (2-12 keV) observations give an 
X-ray luminosity of L I =6.8xl0 37 ergs -1 , and the i?-band 
brightness was ~ 1 mag brighter than its quiescent value. 
The BATSE (20-350 keV) count rate was at least a factor 
of 4 higher than in April/May 1995. Shahbaz et al. (1999) 
determined the true radial velocity of the secondary star 
(#2=215. 5±2.4 kms -1 ) in 1998 May/June, when the source 
was finally in optical quiescence. The only X-ray quiescent 
observations were obtained during March 1996 using ASCA 
(1-10 keV; Robinson et al., 1997) which gave L I =2xl0 32 
ergs -1 . 

In Figure 3 we show the observed radial velocity am- 
plitudes relative to the quiescent value as a function of the 
observed X-ray luminosity at the time of the measurements. 
We have converted the X-ray luminosities, which were ob- 
served with different instruments, into a common energy 
range (0.4-10 keV) using a hydrogen column density of 
Nh = 0.89 x 10 22 cm -2 and a photon power-law model with 
indices 2.8 and 1.5 for the X-ray high and quiescent states 
respectively (see Table 3; Zhang et al., 1997; Robinson et 
al., 1997; Hameury et al., 1997). This energy range is where 



we expect the total radiated power for X-ray transients in 
both outburst and quiescence to lie (Chen, Shrader & Livio 
1997). Note that there is a correlation between X-ray lu- 
minosity and the observed radial velocity semi-amplitude; 
the higher the X-ray luminosity the larger the observed ra- 
dial velocity semi-amplitude, exactly as expected. We can 
use our model to estimate the X-ray luminosity at the time 
when Orosz & Bailyn (1997) took their radial velocity mea- 
surements. We find L x ~ 5 x 10 35 ergs -1 which is consistent 
with the BATSE upper limit. 



6 IRRADIATION OF THE SECONDARY STAR 

It has been known for some time, especially in studies of 
dwarf novae and polars, that substantial heating of the 
secondary star shifts the effective centre of the secondary, 
weighted by the strength of the absorption lines, from the 
centre of mass of the star. This results in a significant dis- 
tortion of the radial velocity curve leading to a spuriously 
high semi-amplitude and a radial velocity curve that may 
be eccentric. Davey & Smith (1992) describe a procedure 
for detecting the effects of irradiation on the radial velocity 
curve of the secondary star, whereby one tests the signifi- 
cance of an eccentricity in the orbital solution. However, it 
should be noted that, although our data does not allow this 
eccentricity test, due to the poor orbital phase coverage, we 
can use the spuriously high radial velocity semi-amplitude 
to show that X-ray heating is present. 

In order to investigate the effects of X-ray heating on 
the secondary star's radial velocity curve we used the model 
described by Phillips, Shahbaz & Podsiadlowski (1999). The 
model uses a crude treatment for X-ray heating, since no 
satisfactory robust model exists for the effects of external 
heating in stars. However, it serves to illustrate the extreme 
effects of X-ray heating. It should be noted that the first or- 
der model of X-ray heating by Brett & Smith (1993), which 
does not include energy transport effects, does show that 
the whole temperature structure of the outer layers of the 
secondary is upset by external heating. Figure 3 shows the 
effects of different amounts of X-ray luminosity on the sec- 
ondary star's radial velocity amplitude. K ohB is computed 
by fitting the predicted curve with an eccentric orbit. The 
regions on the secondary star that are heated do not con- 
tribute to the absorption line flux. The maximum possible 
change that irradiation can have on K b s , based purely on 
geometry, is 15 per cent. However, from our data presented 
in this paper, we observe A#2/#2=0.30±0.05, which when 
compared with maximum possible value based on geometry, 
is significant at the 3-a level. 

In Figure 4 we show how much of the secondary star's 
surface needs to be heated in order to produce the observed 
radial velocity amplitude. We find that based purely on ge- 
ometry 35 per cent of the secondary star's surface is di- 
rectly heated by X-rays produced at the compact object. 
(This fraction only depends on the shape and size of the 
secondary star, which in turn is determind by the q. Us- 
ing the extreme values for q (Shahhaz et al., 1999), we find 
that this fraction changes by less than 1 per cent.) How- 
ever, in order to produce the observed large radial velocity 
semi-amplitude, 60l® per cent of the secondary star needs 
to be heated. The l-o uncertanties quoted here were esti- 
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Table 3. Radial velocity and X-ray luminosity measurements 



K obs L x power-law index Corrected L x 

(kms- 1 ) (ergs" 1 ) N h = 0.89 X 10 22 cm" 1 (erg s" 1 ) 

215.5±2.4 2.0 x 10 32 ASCA (1-10 kcV) 1.5 2.1 x 10 32 

232±2 < 2.4 x 10 36 BATSE (20-350 kev) 2.8 < 1.6 x 10 37 

279±10 6.8 x 10 37 RXTE (2-12 keV) 2.8 8.5 x 10 37 



nated using the 1-a uncertanties in AK2/K2- This result 
may seem surprising at first, since one expects only the re- 
gions of the secondary star facing the compact object to be 
irradiated and yet our result implies that some of the regions 
not directly seen by the compact object are also affected by 
irradiation. However, one should note that effects such as X- 
ray scattering and irradiation-induced flows on the surface 
of the secondary star (Phillips & Podsiadlowski 1999) can 
increase the fraction of the secondary star that responds to 
the X-ray source. Note that the regions on the secondary 
star that are shadowed by the accretion disc will be indi- 
rectly heated by such mechanisms. Therefore K a h s can be 
larger than that expected from heating the inner face of the 
secondary star alone. 



7 DISCUSSION 

The existence of circulation in rotating stars was first proved 
in 1924 by von Zeipel (von Zeipel 1924). He demonstrated 
that for a rotating homogeneous star, the radiative trans- 
port equation and equation of conservation of energy cannot 
be fulfilled simultaneously. This results in the formation of 
meridional motions. In order to maintain a stationary state 
as assumed, one has to demand that these meridional mo- 
tions contribute to the energy transport. In the case of an 
irradiated rotating star, the situation is still more compli- 
cated, since the radiation will induce additional circulation 
currents. 

Evidence for the existence of significant irradiation- 
driven circulation is provided by several sources. For exam- 
ple, the analysis of the optical light curve of HZ Herculis has 
shown this to be heated by its accompanying X-ray source 
HER X-l. Although the main features of the optical light 
variation are well understood (HZ Her is bright when the X- 
ray source is in front of it, its brightness is reduced during the 
occultation of the X-ray source by the secondary), the min- 
imum at phase 0.0 is sharper than expected and indicates 
some additional source of optical radiation at this phase. 
Strittmatter et al. (1973) tried to explain this via the illu- 
mination of the disk by HZ Her. Other attempts were made 
by Pringle (1973) and Bahcall, Joss & Avni (1974). How- 
ever, the most successful explanation was due to Kippen- 
hahn & Thomas (1979). They estimated the energy trans- 
ported from the X-ray illuminated part of the stellar surface 
to the shadowed side, and demonstrated that the minimum 
at phase 0.0 could be reasonably well accounted for (X-ray 
heating without horizontal transport leads to a flat mini- 
mum at phase 0.0). 

In addition, Schandl et al. (1997) found circulation to 
be necessary in order to accurately model the optical light 
curve of CAL 87, an eclipsing supersoft X-ray source. They 
calculated the light curve based on the assumption that an 



accreting, steadily burning white dwarf irradiates the accre- 
tion disk and the secondary star, as suggested by van den 
Heuvel et al. (1992). A simple description of energy trans- 
port on the secondary surface was used and then integrated 
over the whole surface, while conserving the total luminosity. 
They found that significant energy transport of the irradi- 
ated flux to non-illuminated parts on the secondary surface 
is required to simulate the observed lightcurve, particularly 
around the primary eclipse, when the shadowed hemisphere 
of the secondary is in view. 

Recent models for irradiation-induced flows in binary 
stars have been computed by Martin & Davey (1995). 
They considered circulation in gently-heated secondary stars 
(where the incident flux is less than the intrinsic flux) . Their 
2-dimensional calculations included the effects of the Cori- 
olis force and showed upwelling of hot material being car- 
ried preferentially towards the direction of rotation of the 
star. They also concluded that all secondary stars should 
show asymmetric heating, because of the presence of Cori- 
olis forces. Phillips (1999) has recently extended the study 
of circulation to 3-dimensions. As well as including the ef- 
fects of X-ray irradiation i.e. the anisotropic heating of the 
irradiated surface, and the effects of surface radiation stress, 
he also considers the large-scale effects of the rotation of the 
system and includes an approximate treatment of the Corio- 
lis force. His results suggest a realistic analysis of the Coriolis 
force is essential for a full description of stellar circulation. 

In order to study the extent of irradiation of the sec- 
ondary star one requires good quality spectrophotometric 
studies throughout an X-ray outburst, during which the level 
of X-ray irradiation and induced heating changes. This will 
allow the surface intensity distribution across the secondary 
star to be mapped (see Rutten & Dhillon 1994 and Davey & 
Smith 1996), from which effects such as irradiation-induced 
circulation or star-spots can be investigated. 
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Figure 1. Top: Doppler average spectrum of Nova Sco 1994 in the rest frame of the secondary star and an F6 star (HR2906). Bottom: 
Nightly averaged spectra of Nova Sco 1994 (first night at the top, last night at the bottom). The spectra have been normalized and 
shifted vertically for clarity. IS indicates the interstellar 6613A line. 
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Orbital phase 



Figure 2. The radial velocity curve of the secondary star in Nova Sco 1994. The dashed curve is a sinusoidal fit to the data points 
marked with a star. The data have been folded on the orbital ephemeris given by van der Hooft et al., (1998) and 1.5 orbital cycles are 
shown. 
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Figure 3. The effects of different amounts of X-ray heating on the secondary star's radial velocity semi-amplitude. We show the fractional 
change in K ohs as a function of L x . K ohs is computed by fitting the predicted curve with an eccentric orbit. The three observed radial 
velocity measurements are also shown. The two lower curves show the effects of X-ray heating based purely on geometry i.e. only those 
elements of area on the secondary star that are directly seen by the X-ray source are irradiated. The top curve (dashed line) show the 
effects of indirect X-ray heating, calculated by extending the radiation horizon as seen by the X-ray source by a further 24 degrees (see 
Figure 4). The effect of this additional heating is to produce a value for A" b s which is 30 per cent larger than expected purely based on 
geometry. 
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x— coordinate (units of binary separation). 



Figure 4. The irradiated secondary star's Roche lobe in the (x — z) plane. The compact object is at coordinates (1,0). We have assumed 
L x = 8.5 X 10 37 erg s — 1 and the extreme geometrical case with a mass ratio of q=0.44, an inclination angle of i=71° and a disc angle 
of 2°. The dense shaded region shows the area (35 per cent) that is irradiated directly by X-rays produced at the compact object; these 
regions do not contribute to the observed absorption line flux. The less dense region shows the area (60 per cent) which must also be 
heated indirectly in order to produce the large observed radial velocity semi-amplitude. 
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